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Variable angle spectroscopic ellipsometry (VASE) has been used to characterize Si x Gei_ z /Ge superlattices (SLs) grown 
on Ge substrates and thick Si x Gei_ x /Ge heterostructures grown on Si substrates. Our VASE analysis yielded the thicknesses 
and alloy compositions of all layers within the optical penetration depth of the surface. In addition, strain effects were 
observed in the VASE results for layers under both compressive and tensile strain. Results for the SL structures were found 
to be in close agreement with high resolution x-ray diffraction measurements made on the same samples. 

The VASE analysis has been upgraded to characterize linearly graded Si x Gei- x buffer layers. The algorithm has been 
used to determine the total thickness of the buffer layer along with the start and end alloy composition by breaking the 
total thickness into many (typically > 20) equal layers. Our ellipsometric results for 1 fxm buffer layers graded in the ranges 
0.7 < x < 1.0 and 0.5 < x < 1.0 are presented, and compare favorably with the nominal values. 


INTRODUCTION 

Characterization of semiconductor superlattice (SL) 
structures is typically done by x-ray diffraction (XRD) 
and transmission electron microscopy (TEM). Both 
XRD and TEM can be used to quantify SL ordering 
and periodicity. In addition, XRD can give the average 
value of the composition. Both techniques have lim- 
itations; TEM is destructive, XRD requires relatively 
thick samples (~1000 A), and neither technique can 
be used to quantify individual interfaces or quantum 
wells. Recently, variable angle spectroscopic ellipsom- 
etry (VASE) has been shown to be a powerful, non- 
destructive technique for the post-deposition character- 
ization of Si x Gei_ x /Si SLs and other multilayer het- 
erostructures [1,2]. In these studies VASE has been used 
to determine layer thicknesses, alloy composition, oxide 
thickness, number of superlattice periods, and sample 
homogeneity. In this current work we will concentrate 
on Ge rich Si x Gei_ x /Ge SLs and heterostructures. 

Graded composition Si x Gei_ x layers are used in 
the base of Si x Gei_ x /Si heterojunction bipolar tran- 
sistors to increase device speed. Graded composition 
Si x Gei_ x layers (graded to achieve x ^ 0.3) are also 
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used as buffers to relieve strain in the growth of n-type 
Si x Gei_ x /Si modulation doped field effect transistors 
(MODFET) structures. In previous VASE analysis of 
MODFETs, only the high energy portion of the ellip- 
sometric spectra was used. Thus, a graded composi- 
tion layer analysis was avoided, since the graded layer 
was buried below the optical penetration depth of the 
probing light at these energies. In this work a simple 
algorithm for the calculation of the Fresnel reflection co- 
efficients of a linearly graded Si x Gei_ x layer has been 
developed. This allowed a graded Si x Gei_ x layer to be 
characterized in terms of its thickness, and Si content 
at the substrate and ambient surfaces. 


BACKGROUND 

To determine the properties of multilayer struc- 
tures the measured ellipsometric angles, (tan^(A), 
cos A(A)}, must be compared to results calculated from 
well defined models. Linear regression analysis (LRA) 
is used to minimize the unbiased estimator, a: 


1 n 

a2 = n _ m _ ! H { ( tan ** ~ tan *?) 2 

i~ 1 


+ (cos Af - cos A^) 2 } (1) 
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Figure 1: E\ (long dash) and E\ -f Ai (short dash) critical point 
energies for relaxed and strained SixGei-* layers. The relaxed 
curves (Ref. 6) are shown as solid lines. The strain induced split- 
ting (dashed lines) was calculated using Eqs. 2 and 3. Fbr x > 0.5 
the substrate is Si and the strain is compressive, and for x < 0.5 
the substrate is Ge and the strain is tensile. 

where n is the number of observed data points and m is 
the number of free parameters used in the model. The 
superscripts e and c refer to the experimentally observed 
data and the corresponding results calculated by the 
model, respectively. In order to construct a model, the 
dielectric function of all the constituent materials must 
be known. The dielectric function of composite lay- 
ers can be determined using effective medium approxi- 
mations (EMAs) [3]. The energy shift algorithm [4] is 
used to interpolate between Si x Ge!_ x dielectric func- 
tions published at discrete alloy compositions [5]. The 
algorithm requires a knowledge of the functional depen- 
dence of the critical point (CP) energies with the alloy 
composition, usually the E Q (x) t Ei(x ), and E 2 (x) CP 
energies. To improve the analysis of Ge rich Si x Gei_ x 
materials the energy shift algorithm has been modified 
to include the effects of the E\ 4- Ai CP energy [6]. 

Fig. 1 shows the E x and E x + A x CP energies for 
strained and relaxed Si x Gei_ x layers. The in-plane 
strain, in the strained Si x Gei_ x layer results from 
the lattice mismatch between the layer and the sub- 
strate. The shifts in the E\ and E\ 4- Ai CP energies 
for a biaxial (001) strain are given by [7]: 

AE, = ±£i + E,,-I(a; + 4EJ) 1/! (2) 

A(£ i+ A,) = ^1 + E / , + I(a; + 4E1) vi (3) 

where the hydrostatic shift, E H , and uniaxial shear, E s , 
are given by: 

E h = 2£ 1 (l -C l2 /C xl )e (4) 

E s = (2/3) 1/2 Df(l + 2C 12 /C n )e . (5) 


Table I: Comparison of target sample structure with that deter- 
mined by HRXRD and VASE. The period is the sum of the Ge 
and SiGe layer thicknesses, and x(avg) is the average silicon con- 
tent in one period. The parameters of the VASE analysis are the 
oxide, Ge and SiGe layer thicknesses, and the silicon content, x, 
of the SiGe layer. 


Sample 

Source 

Period 

(A) 

i(avg) 

(%) 

d{ Ge) 

(A) 

d(SiGe) 

(A) 

X 

(%) 

HA57 

target 

178 

6.2* 

128 

50 

22.0 

HA57 

HRXRD 

201.1 

8.0 

141. 8 b 

59. 3 b 

26. 6 b 

HA57 

VASE C 

202.8 

8.4* 

127.5 

75.3 

22.7 





±3.9 

±3.5 

±1.6 

HA58 

target 

201 

7.8* 

142 

59 

26.6 

HA58 

HRXRD 

202.1 

8.3 

141. 8 b 

60.3 b 

27. 5 b 

HA58 

VASE d 

204.5 

8.4* 

126.9 

77.6 

22.1 





±4.7 

±4.4 

±2.0 


^Calculated from d(Ge), d(SiGe), and x. 

Estimated from the period, x(avg), and shutter opening/ 
closing times. 


C cr — 0.0137, d(oxide) = 17.9 A ± 0.4 A. 
d <7 = 0.0167, d(oxide) = 25.5 A ± 0.6 A. 


Cij are the elastic stiffness constants, E\ is the hydro- 
static deformation potential, and Df is an intraband de- 
formation potential for the A 3 valence band for a [001] 
uniaxial strain. The energy shift algorithm assumes the 
functional dependence of the relaxed E\ and E\ 4- Ai 
CP energies shown in Fig. 1. Therefore, when modeling 
strained Si x Gej_ x layers the strain induced shifts in the 
E\ and E\ 4- Ai will cause the value of the Si content, 
x y to be over- or under-estimated for a compressive or 
tensile strain, respectively. 

SiGe/Ge SUPERLATTICE STRUCTURES 
ON Ge SUBSTRATES 

Two fifty period Si x Gei_ x /Ge SL structures, HA57 
and HA58, were grown by Si molecular beam epitaxy 
(MBE) on Ge (100) substrates at Hughes Research Lab- 
oratories. The samples were first characterized by high 
resolution x-ray diffraction (HRXRD) [8], the results 
are shown in Table I along with the target structures. 
Each sample was then measured by VASE at several 
angles of incidence (69°, 73°, and 77°) which were cho- 
sen to increase the sensitivity of the ellipsometric angles 
to the sample structure parameters [9]. The measured 
{tan\k, cos A} spectra for sample HA57 are shown in 
Fig. 2 along with the spectra generated from the best fit 
model. Due to the strong absorption in the low wave- 
length region, discrepancies between the measured and 
calculated spectra in this region are a result of devia- 
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Figure 2: Comparison of measured VASE data (symbols) with 
that determined from the best fit model (lines) for the S\ x Gei- x / 
Ge SL sample HA57. Sample HA58 shows similar fitting. The 
structure parameters found from the best fit model are shown in 
Table I. 


tions in the model that occur near the surface. Simi- 
larly, discrepancies in the fitting in the long wavelength 
region are most likely due to problems in the model 
which occur deeper in the structure. In this ease the 
model assumes that each period of the SL is identical. 
Some slight oscillations can be seen in the long wave- 
length region of the cos A spectra which are not repro- 
duced by the model. This indicates that the periods 
of the sample are not identical and that the thickness 
and/or x(avg) of each period may be fluctuating about 
some average value. 

The agreement between the HRXRD and VASE re- 
sults is quite good, especially when comparing the pe- 
riod and x(avg). There is, however, a substantial dif- 
ference between the individual layer thicknesses and x. 
The SijGei-x layer is under a tensile strain as a result of 
being sandwiched between the thicker Ge layers. This 
tensile strain accounts for the lower x value obtained 
by VASE. The differences in the layer thicknesses may 
be a result of interfacial layers which are not currently 
accounted for in the model. 


Nominal VASE 



Oxide 

Si x Gei- x (49%) 
GeO (51% ± 11 %) 

36.5 A 
±6.5 A 

200 A 

S*0.3 7 

SixGe,.* 
x = 29.9% ± 0.8% 

206.1 A 
± 10.1 A 

40 A 

Ge 

Ge 

34.1 A 
± 10.7 A 

1 pm 

Stepped 

SiGe 

Buffer 

Si Substrate 

SiyGei. y Substrate 
y = 25.1% ± 5.2% 
— a = 0.0284 



Figure 3: Structure of sample CL141 showing the target structure 
parameters and those determined from the VASE analysis. 


SiGe/Ge HETEROSTRUCTURES 
ON Si SUBSTRATES 

Two Si x Gei_ x /Ge heterostructures, CL141 and 
CL 171, were grown on Si substrates at the University 
of California at Los Angeles. The nominal structures 
for sample CL141 and CL171 are shown in Fig. 3 and 
Fig. 4, respectively. The Si content of the stepped 
buffer was changed from 100% to 20% or 30% in ap- 
proximate steps of 25%. Therefore, it would require six 
parameters to completely characterize the buffer layer. 
Characterization of the buffer layer is further compli- 
cated by the fact that it is buried deep in the structure. 
As a result sample CL141 was modeled over the range 
300 nm to 550 nm to ensure that the top step of the 
buffer would act as the substrate. Over this range the 
maximum penetration depth is -1000 A and -200 A 
in Sio.3Geo.7 and Ge, respectively. The oxide has been 
modeled as a mixture of Si x Gei_ x ( x same as in under- 
lying layer) and Ge0 2 using the Bruggeman EMA; this 
mixture simulates a surface roughness. The surface of 
these samples is expected to be rough due to the large 
lattice mismatch between the Ge layer and the under- 
lying Si substrate (—4%). The VASE results, shown in 
Fig. 3, agree with in the 90% confidence limits with the 
nominal structure. 

The second sample, CL171, was modeled over the 
spectral range 300 nm to 760 nm. Over this range 
the Si x Gei_ x layer will act as the substrate due to the 
increased thicknesses of both the Ge and SiGe layers 
The measured VASE data is shown in Fig. 5 along with 
spectra generated from the best fit model; the resulting 
structure parameters are shown in Fig. 4. Modeling the 
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Nominal VASE 



Oxide 

Ge (38%)/ 

Ge02 (62% ± 9%) 

36.0 A 
± 2.6 A 

1000 A 

Ge 

S i y Gq _y 
y = 5.8% ± 0.3% 

1263 A 
± li A 

1200 A 

Si 0.2 Ge 0.8 

Si x Gei_ x Substrate 
x = 16.5% ± 1.5% 


1 pm 

Stepped 

SiGe 

Buffer 

a = 0.0269 



Si Substrate 


Table II: Results of VASE analysis of graded Si x Gei_ x layers on 
Si substrates. The nomimal thickness of all samples Is I /im. The 
samples were linearly graded from 100% Si at the substrate to 
x n (nom) at the surface. The parameters of the VASE analysis 
are the oxide thickness, buffer layer thickness (£>), and the Si 
content at the surface (x n )- The SI content at the substrate was 
held constant at 100%. A value of n = 30 was used to model all 
three samples. 


Sample 

a 

d (oxide) 

(A) 

D 

(/im) 

Xn 

(%) 

in(nom) 

(%) 

HA82 

0.0522 

44.8 ± 0.9 

1.04 ±0.04 

76.7 ±0.6 

70 

HA83 

0.0524 

59.3 ±1.1 

1.02 ±0.03 

66.7 ±0.6 

50 

HA17 

0.0202 

43.0 ±0.4 

0.97 ±0.01 

58.7 ±0.3 

50 


Figure 4: Structure of sample CL171 showing the target structure 
parameters and those determined from the VASE analysis. 



Wavelength (rim) 


Figure 5: Comparison of measured VASE data (symbols) with 
that generated from the best fit model (lines) for the thick 
Si*Gei_ x /Ge heteroetructure sample CL 171. The best fit model 
is shown in Fig. 4 along with the resulting values of the structure 
parameters. 


Ge layer as Si x Gei_ x resulted in a ~25% decrease in 
a. Since the Si content of this layer should in fact be 
0%, the value of 5.8% can be attributed to a compres- 
sive strain which results in an over-estimate of the Si 
content. 


SiGe GRADED LAYERS 

Three continuously, linearly graded Si x Gei_ x layers 
were grown by Si MBE on Si (100) substrates at Hughes 
Research Laboratories; the nominal structures are de- 
scribed in Table II. Graded layers are simulated in the 
VASE model by breaking the layer in to n± 1 sub layers. 
The thickness of the i-th sub layer is: 


( D/2n i = 0, n 
\ D/n 1 < i < n 


( 6 ) 


where D is the total layer thickness. For a linearly 
graded layer the Si content of the i-th layer is given by: 


(n - i)xo + ix n 

x t = (7) 

where xo and x n are the Si content at the substrate 
and surface, respectively. As n —► oo this model will 
more closely approximate a continuously graded layer. 
However, values of n > 20 were found to yield nearly 
identical results in the VASE analysis for the graded 
samples used in this study. 

The graded samples were measured by VASE at 
three angles of incidence: 70°, 75°, and 77°. Fig. 6 
shows the experimental data for sample HA17 along 
with the ellipsometric angles generated from the best 
fit model. The fitting is nearly perfect in the cos A 
speotrum and in the long wavelength region of the tan ^ 
spectrum. The poor fitting in the low wavelength region 
of the tan 'F spectrum is most likely due to poor surface 
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Figure 6: Comparison of measured {tan'!', cos A} spectra (sym- 
bols) for the graded Si z Gei_ z sample HA 17 with that determined 
from the best fit model (lines). The structure parameters deter- 
mined from the model are shown in Table II. 

quality. Because the graded layers are uncapped, the 
surface region of the graded layers are expected to be 
strained which can cause the surface to become rough. 
Table II summarizes the results of the VASE analysis. 
It can be seen from this table that the thickness values 
are very close to the nominal values. The values for 
the Si content at the surface (x n ), however, are con- 
sistently much higher than the nominal values. These 
discrepancies can be caused by the compressive strain 
in the surface region of the layer which would result in 
an overestimate of x. The growth log notes that the 
surface morphology of sample HA83 is poor when com- 
pared with sample HA82. This may be the cause of the 
larger discrepancy in x n for sample HA83. 

CONCLUSIONS 

VASE has been used to characterize Si x Ge\- x /Ge 
SL structures grown on Ge substrates. The values of 
the period and x(avg) determined by VASE show excel- 
lent agreement with those same values determined by 
HRXRD. Due to strain effects or interfacial layers not 
accounted for in the VASE models, the individual layer 


thicknesses and x determined by the two methods do 
not agree as closely. Thick Si x Ge 1 _ x /Ge heterostruc- 
tures grown on Si substrates have also been character- 
ized by VASE. Results show close agreement with the 
nominal structure. In addition VASE analysis provides 
a qualitative assessment of the surface roughness and 
strain in the various layers of the structures. Finally, 
VASE has been used to characterize continuously, lin- 
early graded Si x Gei_ x layers in terms of the layer thick- 
ness and Si content at the surfaces. Results agree ex- 
tremely well with the target structure, especially when 
strain in the surface region of these uncaped graded lay- 
ers is taken in to account. 
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